CERN-TH/97-259 



Supersymmetric Q-balls as dark matter 



Alexander KusenkoQ and Mikhail ShaposhnikovQ 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



Abstract 



Supersymmetric extensions of the standard model generically contain stable non-topological 
solitons, Q-balls, which carry baryon or lepton number. We show that large Q-balls can 
be copiously produced in the early universe, can survive until the present time, and can 
contribute to dark matter. 
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Introduction 



The nature of the missing matter in the universe remains one of the most intriguing out- 
standing problems in particle physics and cosmology. Stable particles hypothesized in the- 
ories beyond the standard model have been considered as candidates for dark matter. In 
the supersymmetric extensions of the standard model, the lightest supersymmetric particle 
(LSP), whose stability may be guaranteed by the conservation of R-parity, can naturally 
contribute to the hidden mass. 

Recently, it has been pointed out that non-topological solitons, Q-balls , are gener- 
ically present in the MSSM g). Q-ball is a coherent state of a complex scalar field, whose 
existence and stability are due to the conservation of some global U(l) quantum number. 
In the MSSM the usual baryon and lepton numbers may play the role of such conserved 
quantity for the Q-balls built of squarks and sleptons, respectively. One can ask, whether 
these objects could have formed in the early universe, and whether their lifetime can be long 
enough for them to survive until present and contribute to dark matter. 

Several mechanisms could have lead to the formation of Q-balls in the early universe. 
Non-topological solitons can be created 0] in the course of a phase transition ("solitoge- 
nesis"), or they can be produced via fusion |^, in a process reminiscent of the big 
bang nucleosynthesis ("solitosynthesis"). Finally, small Q-balls can be pair-produced 
at high temperature. None of these scenarios, however, seems to be capable of producing 
large enough Q-balls to survive the subsequent evaporation |^ in the presence of massless 
fermions that carry the same global charge. One might think, therefore, that, in the absence 
of a viable mechanism that could naturally produce very large solitons, no relic Q-balls could 
remain in the present universe. 

In this paper, we show that there exists an effective mechanism for the production of 
Q-balls with enormous charges, big enough to survive until present. The basic idea is the 
following. A complex scalar field (p inside the Q-ball is in a coherent state with time- 
dependent phase: if = i?(x)e"^*, where the radial component -R(r), r = v^^, tends to zero 
as r — s> oo. An infinite size Q-ball with R{r) = const, Q-matter 0, is similar in nature to 
the coherent states of a scalar condensate often encountered in inflationary cosmology. For 
example, formation of the scalar condensate with non-zero baryon number is the starting 



1 



point for the Affleck-Dine scenario for baryogenesis |]T0|- In this scenario, a combination of 
squarks and sleptons, or some other fields carrying a baryon or lepton number, has a large 
expectation value along some flat direction of the potential at the end of inflation. At large 
VEV, the baryon number can be strongly violated by the high-scale physics. As a result 
of the baryon number non-conservation, along with the CP violation, the scalar condensate 
acquires a baryon number. The subsequent evolution leads it into the domain of conserved 
baryon number. From this point on, it can be thought of as Q-matter. 

Depending on the dynamics, an initially spatially-homogeneous scalar condensate that 
carries a conserved U(l) charge may become unstable with respect to small coordinate- 
dependent perturbations and develop a spatial pattern that comprises domains of high and 
low charge density. If a theory admits non-topological solitons fl], the lowest energy state 
in a given charge sector is a Q-ball. Very large (and, thus, stable) Q-balls can, in fact, be 
produced this way. 

A related issue is that of the string moduli. If their evolution leads to formation of 
isolated solitons, it can help alleviate the cosmological moduli problem |11, l^- The weakly- 
interacting moduli may possess some energy density during or after the nucleosynthesis. 
However, if this energy density is incarcerated inside Q-balls, it may be harmless. 

In what follows, we first discuss the stability and decay of Q-balls. Then we analyze 
formation of Q-balls through solitosynthesis and collisions. We then consider the evolution 
of a complex scalar field and describe the formation of a spatial pattern that can evolve into 
Q-balls. Finally, we discuss the contribution of relic Q-balls to dark matter and comment 
on the possibility of solving the cosmological moduli problem through the incarceration of 
malicious moduli. 



1 Lifetime of Q-balls 

The fate of Q-balls in the early universe is ultimately determined by their lifetime. If the 
soliton in question is made of a scalar that has no interactions with light fermions, then 
it is entirely stable. If such objects had formed through the breakdown of a homogeneous 
condensate, they would have survived until the present and would contribute to the matter 
density of the universe as a form of dark matter. 
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If, however, the scalar particles can decay into lighter fermions that carry the same global 
charges, as is the case with the squarks and sleptons of the MSSM, then it is the evaporation 
rate that determines the lifetime of a Q-ball 0. Clearly, it is only the fermionic decay 
modes that ought to be considered in this respect, because, in the sector of scalar degrees of 
freedom, the Q-ball is the state of minimal energy, and it cannot decay into bosons. 

Depending on the scalar potential [/(yji, the mass of a large Q-ball may grow 
with charge as ~ Q^, where < p < 1. If the quantity 

where qi is the charge of the field (pi (some of the charges may be zero), is minimized at a 
finite value of yj, then the mass of a large Q-ball is proportional to the first power oi Q [Q. 
In this case the mass per unit charge does not depend on the size of the Q-ball. 

However, if U {(p) is essentially flat, i. e. if it grows slower than the second power of then, 
in general, p < 1 |13[. Such potentials can arise naturally in theories with supersymmetry 
breaking communicated at low energy, for example, by gauge interaction (for review, see, 
e. g., Ref. and references therein). If the potential is essentially flat for large (f ^ m^, 
then ~ Q^^^ [|13|. We recapitulate briefly the argument of Ref. |jl3|. The field inside the 
Q-ball is ipi{x,t) = ipi{x) exp{iqiUjt} |^. The value of u is determined by extremizing the 
functional 
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ujQ, (2) 



where U^{ip) = U{ip) — ^u;^ J2k Ik Iv'fcP- Written in terms of dimensionless variables ^ = ujx 
and ipi = (fi/uj, ^ auj + h/uo^ + uoQ, if U{(p) ~ const. Here a and h are independent of uo. 
Therefore, the extremum of functional @) is achieved for uj oc Q^^^, which corresponds to 

oc (3) 

By virtue of the relation (|^), a sufficiently large Q-ball is entirely stable against decay 
into all but massless fermions because its energy per unit charge can be less than the mass 
of the lightest fermion. In the MSSM, there are Q-balls (B-balls) made of squarks and Higgs 
fields [^, for which the baryon number plays the role of the U(l) charge discussed above. 
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The lightest baryons, nucleons, have mass of order m„ ~ 1 GeV. A B-ball is entirely stable 
if it has charge 



Q = B>i-^ ~ 10«, (4) 
Vm„/ 

where we took ~ 10^ GeV for the squark mass. 

If the Q-ball has a lepton charge, but no baryon charge, then it can evaporate |^ by 
emitting lepton number from the surface with (nearly massless) neutrinos. It is beyond the 
scope of the present work to evaluate the decay width of slepton Q-balls in the MSSM taking 
into account the variety of fermionic decay modes. Instead, we will obtain a rough estimate 
by adopting the discussion of Ref. p to our case. 

The model of Ref. |p had a Yukawa coupling of a scalar field (f to massless fermions. 
Inside the Q-ball, the Dirac sea of fermions is filled; and the Fermi pressure inhibits further 
decay. The decay of the scalar condensate can, therefore, proceed through the surface only, 
and its rate is suppressed by the surface-to-volume ratio for large Q-balls. The rate of such 
decay is 



dA dt 1927r2 ' 

where A is a Q-ball surface area. For the lifetime of an L-ball to exceed the age of the 
Universe, its charge should be Q > 10'^^. 

The effects of gravity on Q-balls are very small and will be neglected, as long as the size 
of the soliton is much larger than the corresponding Schwartzschild radius. For a Q-ball in 
a flat potential, this is the case ii Q <^Qg = {nip/m^Y ~ 10^^ ~ 100 GeV. A Q-ball 

with charge that exceeds Qg would be a black hole. It would evaporate through Hawking 
radiation, and its global charge would be lost. In our analyses we never encounter Q-balls 
with such large charges, and we don't know of any mechanism by means of which they could 
form. 

If the potential is not fiat, then ~ Q 0, and the Q-ball would decay through 
evaporation [P] regardless of the value of its charge. (The charge needed to survive the 
evaporation would exceed the black hole limit by many orders of magnitude.) 

At finite temperature, the existence of Q-balls depends on the form of the temperature- 
dependent effective potential. If it admits Q-balls, then the thermal processes do not affect an 



existing soliton. If, however, a Q-ball is formed out of equilibrium at some temperature T in 
a system that does not allow Q-ball solutions, the collisions with particles in the hot plasma 
can lead to an erosion of the condensate. For large solitons this process is strongly suppressed 
by the small penetration depth of the outside particles inside the Q-ball. A particle ip that 
interacts with the field ip receives a contribution to its mass of order m^(x) ~ g{(f{x)) inside 
the Q-ball; outside Q-ball ~ gT. At temperature T a. ip particle cannot penetrate inside 
the Q-ball beyond the point x where f{x) ~ T/g. The condensate can be thought of as a 
coherent state of particles with mass rricond ~ rn^/Q^^^ and density ujLp^{x). If a Q-ball is 
created at the time tj, the total number of collisions with the particles in plasma during 
the lifetime of the cooling universe is 



where a ~ g"^ /T"^ is a cross-section of particle-particle collision, the velocity of the ip particle 
at a given point is determined by energy conservation. The stopping point, tq is defined by 
the relation m^{r) = T, or g(f{x) = T. For squark B-balls (7 ~ 1. One can use = Mg/t 
to relate temperature to time in a 
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radiation dominated universe. Here and below Mq = Mp / {1.66 njjrj-), and Uefj is the 
number of effective degrees of freedom. In each collision, only a fraction of energy, 6E ~ 
Tirricond/f^^Y , is transmitted to the condensate. Therefore, using (|^), the total energy 
transmitted to the Q-ball is 



This is sufficient to knock out AQ ~ AE/m^p particles from the condensate. For a "flat" 
potential. 



2 Solitosynthesis and collisions 

Solitosynthesis of non-topological solitons |], is the process of charge accretion that 
may take place in a system at finite temperature with non-zero charge asymmetry rj. Q-ball 




(6) 
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A Q-ball can survive the erosion by thermal plasma if AQ < Q. 
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is the minimum of energy E in the sector of fixed charge, but is also a state of small entropy 
S (because it is an extended object). For high temperature it may not be the minimum 
of free energy F = E — TS because of the TS term. At some small enough temperature 
To ~ In^^l, however, the second term in F is no longer important, and the gain from 
minimizing E overwhelms the loss from lowering the entropy §|, 0- At that point, Q-ball 
is the minimum of free energy. In the absence of light charged fermions, a copious production 
of Q-balls can take place. 

However, if there are light fermions carrying the same global charge as the scalars, then 
the charge asymmetry in the state of the lowest free energy is accommodated entirely in the 
fermionic sector, thus making solitosynthesis impossible. For example, the solitosynthesis of 
baryonic and leptonic Q-balls in the MSSM is only possible in a B or L breaking minimum 
of the potential, where the quarks or leptons are sufficiently massive 0. 

Even in the absence of light charged fermions, solitosynthesis is unlikely to produce large 
Q-balls unless the charge asymmetry is very large |^. The baryon asymmetry of the 
universe is as small as 10"^'', which is certainly not sufficient for synthesizing large enough 
B-balls to satisfy the stability bound @. One can ask whether the subsequent mergers of 
small Q-balls can increase the average size considerably. 

Regardless of the mechanism that lead to the formation of Q-balls, their size can further 
change because of the collisions in which two Q-balls would merge and form a larger soliton. 
For simplicity we will neglect the evaporation of Q-balls assuming that they are made of 
stable scalars. The effect of collisions is described by the kinetic equation 



(9) 



where A''^ is the density of Q-balls with charge Q, and H is the Hubble constant. We note in 
passing that if the individual particles are included as Q = 1 objects, the same equation (|) 
can describe the charge accretion via absorption of particles by Q-balls. The cross-section 
for the merger of two Q-balls is determined mainly by their geometrical size, ~ 
and the average velocity of their non-relativistic Brownian motion is f ~ ^JT /m^. For a 
thick-wall Q-ball in a flat potential the quantity {(J^v) oc -R^/ .Jfn^ changes with charge very 
slowly, as Q^^^. Therefore, for a large range of charges, it can be considered constant. One 
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can define the reduced number density of Q-balls as = N^/T^. Tlien 

d 



J ngnQ^qdq - 2nQ J Uqdq, (10) 



wliere r = Tmax[^ — {T/Tq^^'^] and varies from r = to T^ax = Mo{'^q)T^^'^ / y/i^- This 
equation can be solved by applying a Laplace transform to disentangle the integrals on 
the right-hand side. The resulting equation is a first-order differential equation for fs{t) = 
J^nQexp{-Qs)dQ, 

^fs = fs-2fJo- (11) 

Let us consider, for example, an initial distribution that peaks near Q = Qq. The value of 
Qq can be set by some earlier process, for example, sohtosynthesis. We require that at r = 0, 
nQ{l) = e6{Q — Qq), and, hence, /^(l) = eexp{— Qo^}. This initial condition corresponds 
to a set of equal size Q-balls with charge Qq each and number density e that start merging 
through collisions at r = 0. The solution is 

The inverse Laplace transform of /s(t) in equation ( |12|) is a contour integral that can be 
found by summing over the residues of the function (/s(t) exp{Qs}) in the complex plane. 



(1 + er) 



er 



1 + er 



{Q-Qo)/Qo 

E^(QM-ri), (13) 



where n is an integer. The distribution (|T^) describes an array of delta-functions at integer 
values of Q/Qq weighted differently for different r. This is to be expected because each 
collision changes the charge of a Q-ball by an integer multiple of Qq. 

As temperature goes to zero, r T^ax, the distribution of Q-balls approaches (for 

^'^max 1) 

lim riQ (t) oc - — ^exp( (14) 

[eTmaxr \ (^^maxQoJ 

for integer Q/Qq. This distribution is suppressed for large values of Q/Qq 3> eTmax- 

This process, which, for reasons explained earlier, is only relevant for weakly interacting 
stable particles, may considerably increase the average size of a Q-ball. 
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3 Instability and pattern formation in the motion of a 
scalar condensate 



We now turn to the mechanism that can, in fact, lead to the formation of very large Q- 
balls in the early universe. In particular, huge Q-balls made of squarks and sleptons can be 
produced this way. 

Evolution of a spatially homogeneous condensate carrying a baryon or lepton number has 
been studied in detail in connection with the Affleck-Dine [|10] mechanism for baryogenesis. 
The scalar under consideration is a complex field that carries some U(l) charge (e. g., the 
baryon number) and can be a combination of squarks and sleptons of the MSSM [|I5[. If the 
U(l) charge is preserved by the low-energy physics, the eventual decay of such condensate 
can contribute to the corresponding charge asymmetry. 

It is possible, however, that the spatially homogeneous condensate develops an instability 
that leads to the formation of domains with higher and lower charge density. If the theory 
admits Q-balls, such pattern can evolve into Q-balls before the condensate is eroded by 
decay. Thermal fluctuations cannot wash out the coherent state as long as the scalar VEV 
is larger than gT, where (7 is a typical coupling. This is because all the particles coupled to 
Lf with a coupling g have masses of order g{ip) ^ T. After the Q-balls form, the subsequent 
decay of the condensate may become impossible (if the Q-ball energy per unit charge is less 
than the mass of the lightest Q-charged fermion) or slow (if the charge is sufflciently large). 
The relic Q-balls can thus survive until present and contribute to the dark matter in the 
universe. 

Under very general conditions the scalar condensate can be described by the classical 
equations of motion in the effective potential with initial conditions determined by the high- 
scale physics or inflation. It is convenient to write the complex fleld = i?e*^ in terms of 
its radial component and a phase, both real functions of the space-time coordinates. We are 
interested in the evolution of the scalar fleld in the low-energy domain, where the baryon 
number violating processes are suppressed, and we will assume that the scalar potential 
preserves the U(l) symmetry: U{ip) = U{R), where U{R) may depend on time explicitly. 
The classical equations of motion in the spherically symmetric metric ds"^ = dt^ — a{tYdr'^ 
are 
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R + 3HR- ^^AR-Q'^R+ ^^(diQ?R+(dU/dR) = 0, (16) 

where dots denote the time derivatives, and the space coordinates are labeled by the Latin 
indices that run from 1 to 3. The Hubble constant H = a/a, where a{t) is the scale factor, 
is equal to t"^/^ or t^/^ for the matter or radiation dominated universe, respectively. 

A spatially homogeneous solution of equations (|T5|) and (|T6|) has vanishing gradients of 
R and Vt. If one neglects the expansion of the universe {H = 0) and if -R = 0, the solution 
describes a universe filled with Q-matter 0. If the theory admits Q-balls, Q-matter can 
become unstable with respect to its break-up into Q-balls. In a static universe, the conditions 
of stability were analyzed in Ref. [0. The results cannot be applied directly to the case 
of the expanding universe, where the equations of motion are time-dependent and the red 
shift of linear perturbations is important. Under some conditions spatial inhomogeneities 
may destabilize a given solution of equations (0) and (|16|). 

From the equations of motion (|15D and (0), one can derive the equations for small 
perturbations 6Q and 6R: 

5n + 3H{5h)-—-Ai5n) + ^{5n) + —{6R)--^5R = 0, (17) 

a [tj ri ri ri 

5R + mi5R)- ^--M5R)-2Rhi5h) + U''5R-tf5R = 0. (18) 
a"'(t) 

In general, the analysis of instabilities must be done numerically. However, one can 
identify the growing modes and estimate the length scale of the resulting pattern by using a 
WKB approximation. This approximation is valid as long as the background solution varies 
slowly with time. 

To examine the stability of a homogeneous solution ^p{x,t) = Lp{t) = i?(t)e*^*^*\ let us 
consider a perturbation 5R,Sil oc and look for growing modes. Re a > 0, where 

a = dS/dt. The value of k is the spectral index in the comoving frame and is red-shifted 
with respect to the physical wavenumber k = k/a{t) in the expanding background. Of course, 
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if the instability develops, the linear approximation is no longer valid. However, we assume 
that the wavelength of the fastest-growing mode sets the scale for the high and low density 
domains that evolve into Q-ball eventually. This assumption can be verified post factum by 
comparison with a numerical solution of the corresponding partial differential equations (|15D 
and (|T^), where both large and small perturbations are taken into account. 
The dispersion relation follows from the equations of motion: 



a + 3Ha + ^ + —a 
R 



^2 

+ 3Ha + — - fi^ + U"{R) 





R 








a = 







If {Cl"^ — U"{R)) > 0, there is a band of growing modes that lies between the two zeros of 
a{k), < k < kmax, where 



k^a.it)=ait),Jn^~U"{R). (20) 

If kmax{t) defined by relation ( ]20|) is constant or growing with time, then each mode has 
an unlimited time to develop and the size of a resulting inhomogeneity in the comoving frame 
is determined by the value of k that corresponds to the maximum of a. 

If, on the other hand, kmax{t) decreases with time, then the modes may not have the time 
to grow sufficiently before they are red-shifted from the resonance. A constant or growing 
with time width of the band (^) is a sufficient (although not necessary) condition for the 
instability. 

Quantitatively, the amplification of a given mode k is characterized by the exponential 

of 



S{k) = / a{k,t)dt. (21) 



We look for the value of the best-amplified mode with wavenumber kf,est that maximizes S{k). 
The instability develops when e^^'''""'*^ ^ 1. The amount of amplification needed in a given 
system is determined by the size of the fiuctuations. These can be thermal, or of a different 
origin (see, e. g., the discussion in Ref. [p!7|). The critical value of S{khest) ~ lia{R/SR) signals 
the non-linear regime and corresponds to the breakdown of the homogeneous solution. 
As an example, let us consider the following potential: 
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f/(^) = m^^ In + cH'^^^ + 4-(^V)', (22) 

\ S / P 

where c and A are some constants and H ~ 1/t. Such potential can arise in theories with 



supersymmetry breaking communicated at the scale m„ <^ nip |18, 19| and is similar to those 



discussed in Ref. |I9[ in connection with baryogenesis. The mass term proportional to 
the Hubble constant appears naturally from the Kahler potential as / (i^^((x^x)/m^) yj^y^ 
1^ , where x is some field with non-zero density, whose F-component has a VEV. This term 



essentially characterizes supersymmetry breaking by the finite energy density that dominates 
the universe at time t. 

We will assume that the universe before Q-balls formation is matter dominated. For large 
the logarithmic term in ( ^2] ) can be neglected and the equations of motion admit an 
adiabatic solution for R(t) ~ and = 17 ~ that follows closely the location of the 

minimum of U. The conservation of charge per comoving volume implies Cl{t)R'^{t)a{tY = 
const. 

It is clear from equation (|20|) that the band of resonance, bounded from above by k^ax ~ 
becomes narrower with time. In other words, each mode in the amplification region is 
red-shifted faster than k^ax, and the amplification may not be sufficient. We note that the 
time dependence of uj is the same as that of H, both are inversely proportional to time. 

For small ~ m^, the first term in equation p^ ) dominates: U{(f) ~ m?^ip'^ip. In this 
case, the equation of motion for R is linear and its solution is ~ 1/t, with a constant 
uo = m^. Clearly, in this regime there is an unlimited growth of all modes in the band of 
instability because k^ax ~ ct(^) is growing with time and the band widens. 

Qualitatively, the behavior of u is shown in Fig. As long as the logarithmic term in 
the potential can be neglected, uj changes as 1/t, roughly, at the rate of the expansion of 
the universe, and its value is of order H. However, at later times uj approaches a plateau 
UJ ~ TTig. A characteristic length of growing instability is determined by the fastest-growing 
mode at the time when the rate of change in uj becomes different from 1/t. 

To estimate a charge of Q-balls produced we assume that Q-matter was formed at some 
time to with charge density go = Q'(^o) = ^^oVo- ^ matter dominated universe g(t) = 
go(to/t)^. Naturally, the wavelength 1/k ~ x/H of the best-amplified mode cannot exceed 
the size of the horizon at the time tj when the instability develops, x < 1- The total charge 



11 



CO 



m 

s 




t 



Figure 1: At early times uj changes at, roughly, the same rate as the expansion of the 
universe. All modes of instability are red-shifted before they can grow sufficiently. However, 
at later times uj approaches a constant value ~ and the instabilities develop. The size of 
a domain is determined by the fastest-growing mode at the time when the rate of change in 
UJ becomes different from 1/t. 

of Q-matter inside the part of horison corresponding to the best amplified model gives an 
estimate for the charge of Q-ball: 



The time ti and parameter x are to be found from numerical calculations. We performed 
them for a potential (|2^) for different sets of parameters. The amplification S{k) of a given 
mode was found through numerical integration. We took the initial time of the Q-matter 
formation to be to = 10"^ GeV-^ For A = 1/2, c = 1, = 10^ GeV and go = lO^^GeV^ 
(this corresponds to (fo = 2.5 x 10^^ GeV and ujq = 1.6 x 10'^ GeV) we observed that the 
best-amplified mode was two orders of magnitude smaller than the horizon size, x — 10^^, at 
the time when the instability developed {Smax{k) — 30 at ti ^ 150 GeV~^). This corresponds 
to the average charge of emerging Q-balls of order 10^^. For a different set of parameters, 

= 10^ GeV, A = 0.5 x 10"^ and go = 10^^ GeV^ (which correspond to ipo = 2.8 x 10" 
GeV and uq = 1.3 x 10^ GeV), we obtain U ~ 1.5 x 10^ GeV~\ x ^ lO^^ and Q ~ 10^°. 
The band of unstable modes is shown in Figure Of course, the size of Q-balls depends 
on a model. We have demonstrated, however, that very large Q-balls can be born from 
an initially homogeneous scalar condensate. In particular, the flatter, than logarithmic, 
potentials produce even larger Q-balls. It would be desirable to study the instabilities beyond 
the linear analyses performed so far, when the solution is signiflcantly different from the 



Q ~ qotlUx^. 



(23) 
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Figure 2: The amplification S{k) = j/™""' a{k,t)dt of growing modes is computed numeri- 
cally for different values of k. The fastest-growing mode has the wave number k ~ lO^H at 
the time when the instability develops. An individual Q-ball gathers charge from, roughly, 
10^^ of the size of the horizon. 

original homogeneous one, as illustrated in Figure 0. We leave this for future work. 

4 Surviving Q-balls as dark matter 

It is beyond the scope of this paper to examine in detail the implications of solitonic dark 
matter formed from the breakdown of a scalar condensate. We will limit our discussion to a 
few general remarks. 

We have seen that B-balls formed through the instability can easily satisfy the bound 
(^). If finite temperature effective action admits B-balls, they will survive till present. If not, 
they may or may not be eroded by the thermal plasma, depending on the model parameters 
and temperature Tj at the time of the Q-ball formation. From equation (H) they are not 
destroyed if Tj < m^p{m^Q / Mq)^ . Surviving B-balls would contribute to the dark matter in 
the Universe. The amount of the dark matter is determined by the efficiency of the conversion 
process that leads from a uniform condensate to B-balls. As the instability develops, part 
of the charge (baryon number) stored in the condensate can be lost in decay and give rise 
to the baryonic matter. It seems reasonable to assume that the amount of baryonic charge 
locked inside the Q-balls can be roughly comparable (within a few orders of magnitude) to 
the amount of charge that leaked out with individual nucleons. 

If the scalar condensate that evolves into Q-balls has a lepton charge, but no baryon 
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Figure 3: The charge density per comoving volume in (1+1) dimensions for a sample 
potential analyzed numerically during the period when the spatially homogeneous condensate 
breaks up into high- and low-density domains. Two domains with high charge density are 
expected to form Q-balls. 

charge, then a Q-ball can evaporate by emitting lepton number from the surface with 
(nearly massless) neutrinos. We have seen that, for the lifetime of an L-ball to exceed the 
age of the Universe, its charge should be Q > 10'^^ according to equation (||). Smaller L-balls 
must have evaporated by now. Decay of the L-balls during or after the nucleosynthesis may 
have important ramifications as it would lead to the entropy increase and, if it happens at 
later times, could distort the spectrum of the cosmic microwave background radiation. This 
may lead to constraints on the parameters of the MSSM, as well as on the initial conditions 
after inflation. 

Our discussion has bearing on the issue of the so called moduli problem |Tl], |T2|. Small 



couplings to matter fields are typical for the moduli fields of string theory. This poses a 
serious problem for cosmology since the energy density carried by the coherent motion of 
the moduli condensates does not decay until after the nucleosynthesis. The products of 
such decay can over-produce the entropy, as well as destroy the nuclei. The mechanism we 
have described may help alleviate the moduli problem by incarcerating moduli inside slowly 
evaporating Q-balls. On the one hand, this slows down the decay, thus allowing for the 
moduli Q-balls to still be present in the Universe and contribute to dark matter. On the 
other hand, whatever decay would have taken place, the detrimental effect of the energetic 
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photons on the nuclei is confined to the close vicinity of Q-balls separated by some vast 
unaffected areas, where the standard nucleosynthesis can take place. 

Let us suppose that in the present Universe there is, roughly, the same (to within a few 
orders of magnitude) baryon number in baryonic Q-balls and in ordinary matter. According 
to our estimates, the charge of an individual stable B-ball may be from 10^ to 10^° or more; 
and the corresponding number density of Q-balls is then 10^^° to 10~^ per nucleon. Perhaps, 
even much larger Q-balls may emerge from the scalar condensate in models different from 
ours. The larger Q-balls from this range may be difficult to detect because they are extremely 
rare. However, the possibility of smaller, and, hence, more populous B-balls as a dark matter 
candidate seems to be very appealing from an experimental point of view. Some of the heavy 
relic Q-balls, attracted by the gravitational fields of stars and planets, can sink to the center 
and remain there. It is conceivable that the deep interior of the small planets might become 
accessible for exploration in the future and reveal storages of primordial Q-balls. Other 
solitons may be present in the interstellar medium. However, since their number density is 
expected to be very small, detection of such objects seems to be a great challenge to the 
experimentalists' ingenuity. It is also of interest to understand the effect of large Q-balls on 
the structure formation in the early universe. 

We thank G. Dvali, H. Murayama and P. Tinyakov for helpful discussions. 
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